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Abstract 

Generalized symmetry energy coefficients of asymmetric nuclear matter are obtained as the 
screening functions. The dependence of the isospin symmetry energy coefficient on the neutron 
proton (n-p) asymmetry may be determined unless by a constant (exponent) Z which depend 
on microscopic properties. The dependence of the generalized symmetry energy coefficients with 
Skyrme forces on the n-p asymmetry and on the density, only from .5 up to 1.5 poi a-re investigated 
in the isospin and scalar channels. The use of Skyrme-type effective forces allows us to obtain 
analytical expressions for these parameters as well as their dependences on the neutron-proton (n- 
p) asymmetry, density and even temperature. Whereas the density dependence of these coefHcients 
obtained with Skyrme forces is not necessarily realistic the dependence on the n-p asymmetry 
exhibit a more consistent behaviour. The isospin symmetry energy coefficient (s.e.c.) increases as 
the n-p asymmetry acquires higher values whereas the isoscalar s.e.c. decreases. Some consequences 
for the Supernovae mechanism are discussed. 
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1 Introduction 



The symmetry energy coefficients are of interest for the understanding of several aspects of nuclear 
structure and also reactions as for instance (neutron rich) nuclei structure, giant resonances and nuclear 
heavy ion collisions at intermediary and high energies |]|, |2|, ^. There are also astrophysical 
motivations for this as, for example, in the supernovae mechanism and r-processes in nucleosynthesis 
[^, 0, ^. However, these coefficients can depend on the medium properties as density, temperature 
and neutron-proton asymmetry being these dependences relevant for different phenomena. Although 
one usually only considers the coefficient of the isospin channel (a,-) there are other coefficients for the 
spin, spin-isospin and scalar channels of the effective nucleon nucleon (NN) interaction. The former 
is usually defined as the energy difference between the unpolarized and polarized nuclear matter (ag-) 
while the latter to the completely symmetric unpolarized and asymmetric polarized neutron matter 
(ao-r)- However these two last parameters are believed to be small and will not be considered in the 
present article. The most well known and studied of these parameters is the isospin one (a,-) and we 
will refer to it unless explicitly quoting the others. This parameter measures the restoring force of the 
nuclear system to a perturbation which separates protons from neutrons. It corresponds to a cost in 
energy which however depends on the asymmetry of the medium. 

In the macroscopic mass formulae the symmetry energy coefficient (s.e.c.) (in the isospin channel) 
measures the difficulty of increasing the neutron-proton (n-p) asymmetry for a stable system. It 
contributes as a coefficient for the squared neutron-proton asymmetry: 

E/A = Ho{A)+ar{N-Zf/A\ (1) 

where Hq does not depend on the asymmetry, Z, N and A are the proton, neutron and mass numbers 
respectively. An analogous form for the scalar channel will be considered yielding a parameter related 
to the nuclear matter incompressibility. Expression (|l]) is obtained from the Fermi gas model Q as 
well as from microscopic calculations ||l^ but it is also known that there are important corrections due 
to nucleon-nucleon interactions and correlations. Higher orders effects of the asymmetry (proportional 
to {N — Z)" for n 7^ 2, n=l, 3) correspond to the explicit breaking of isospin symmetry among other 
effects |pl|. They are expected, in principle, to be less important for the equation of state (EOS) of 



nuclear matter |12, 13|. However we show that there may occur relevant deviations for n-p asymmetries 
by studying the generalized s.e.c. as the nuclear matter screening functions. 
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For the purpose of comparison we quote one work in which the density dependence of the isospin 



s.e.c. worked out |14]. In this reference the relevance of for the equation of state (EOS) of dense 
neutron stars was studied with an effective form for the nuclear interaction. In that work the following 
parametrization was obtained: 

ar{p) = S{p) = (2i - l)^E%{ul - F{u)) + 5oF(n), (2) 

where u = p/ Pq, po being the saturation density of nuclear matter, Ep is the Fermi energy of nuclear 
matter, 5*0 is the Ur at the saturation density (nearly 30 MeV) and F{u) is a generic function which 
satisfies the condition F{1) = 1. With non relativistic microscopic approaches, a-^ is found to vary 
nearly linearly with the density until p ~ 2.5po ^-nd to saturate at higher densities. However, this 
behaviour may depend on the nucleon-nucleon potential |1C, 15, |l^ ]. 

In the present article we perform a detailed investigation of the nuclear matter s.e.c. as the 
screening functions of nuclear matter for the isospin and scalar channels. Preliminary results were 
presented in Now we extend that calculation for the case of Skyrme forces with more involved 
density dependence [|l8|, |l9| and make a quite complete analysis of them. Different effective Skyrme 
forces are considered in order to assess the possible behaviour of these functions. The use of these forces 
allows us to derive analytical expressions for the generalized s.e.c. One may argue that Skyrme effective 
interactions are parameterized purely on phenomenological grounds and therefore would not have 
predictive power for nuclear matter at densities different from po (the saturation density) or at very 
large n-p asymmetries. In the present work, however, we consider different Skyrme forces fitted from 
diverse sources including one derived from studies with realistic interactions | |2[)[ |. Therefore, a study 
with these different forces can supply some information on the behaviour of such phenomenological 
models. Furthermore, as Skyrme parameterizations are frequently used, this kind of information can 
be relevant in diverse situations. Besides that, a general calculation is done which yields a general 
expression for the isospin dependence of the screening function. The dependences of these generalized 
coefficients on the n-p asymmetry and on the nuclear density (in a short range: from O.Spo up to 1.5po) 
will be studied in this article. The relevance of some results for the Supernovae (SN) mechanism is 
discussed. In the next section the general expression for the screening functions is derived and the 
particular expression for Skyrme forces is presented. The specific form of the functions of the Skyrme 
parameters is shown in the appendix for each of the channels. In section 3 the results are presented 
and discussed. In the final part the conclusions are summarized. 
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2 Generalized Symmetry Energy Coefficients 

Generalized nuclear matter symmetry energy coefficients will be investigated in the following. It is 



interesting to review and to extend a qualitative argument from |21| for exploring them. 



2.1 General Remarks 

Let us consider a small amplitude (e) external perturbation which acts, through the third Pauli 
isospin matrix t^, in nuclear matter separating nucleons with isospin up and down This originates 
a fluctuation 5p = pn — Pp of the nucleon density. The total energy of the system can be written as: 

H = Ho + ar— + e6p, (3) 
P 

where Ot is the isospin symmetry coefficient and Hq does not depend on the asymmetry. In the 
equilibrium the following stability conditions must hold: 

2^ + 6 = 0. 4 
dp p 

The ratio of the amplitude of the fluctuation generated by the external perturbation to the amplitude 
e of this perturbation yields the static polar izability (11) (now generalizing for any channel as done in 



22, ^3] with (s,t) for spin, isospin numbers): 



^ - n-* = -J^, (5) 

Where the Ag^t are equal to the s.e.c. In particular in the isovector channel: ^o,i = o-t- This expression 
corresponds to the static limit of the response function of symmetric nuclear matter | ]21| |. 

Now let us consider asymmetric nuclear matter. In the usual parameterizations for the mass 
formulae, there are several effects which depend on the induced asymmetry 5p = pn — Pp = /?• 
One can thus extend expression (^) to include a more general dependence on the asymmetry. Our 
approach consists in considering a general parameter A which is also a function of the n-p asymmetry 
(and therefore a function of /?) instead of the s.e.c. as usually defined in mass formulae of the type of 
(^). If the same external perturbation is introduced in the energy density we can write: 

= i7o + ^°''(/3)— + e/3. (6) 
P 

^ This argument is valid for all the four channels -isovector, spin, spin-isospin and scalar- with suitable modification. 
It is enough to consider other external perturbations: a, crara and 1 for the spin, spin-isospin and scalar channels 
respectively. 
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It is worth to note that A^'^ (in a notation of a generic channel) has rather a dependence on the n-p 
asymmetry b = pn/ Pp — ^- The exact form of this dependence is not the same that the dependence on 
the fluctuation In other words we could write, for the sake of completeness, A = A{(3, b). However 



these two parameters should be related to each other. In [23| two different prescriptions were discussed 
for P in the calculation of the response function of asymmetric nuclear matter. The most reasonable 
one leads to the following relationship between the fluctuation /? and the asymmetry b: 



/5 = spn {m) ' 



Where 6pn is the neutron density fluctuation. In the symmetric limit /3 = 2pn and in the opposite limit, 
in neutron matter, /? = 5pn, as expected. Considering this, we are lead to expect that the dependence 
of A on the asymmetry b is related to its dependence on the density fluctuation /3. This allows us to 
write shortly: A = A{P). It is important to stress that the prescription (^) is model-dependent and 
different choices for it yield different results for the static screening functions as shown below. The 
dynamic response functions are less sensitive to this prescription, but not completely independent [p3| . 
The above prescription (^) is based on the assumption that the density fluctuations are proportional 
to the respective density of neutrons and protons, i.e., Spn/f3 = Pn/p, being p the total density. In 
spite of being rather well suited for the isovector channel, this kind of assumption can be considered as 
a simplified model for the other spin/isospin channels in asymmetric nuclear matter. The derivative 
of A with relation to /3 is therefore related to dA/db. 

If the system is stable with relation to the density fluctuations induced by the external source, the 
following equation holds: 

§.2^/?^|£i?! + , = „. (8) 

op p op p 

By substituting the variables 11 = /3/e and B = A/ p va. this equation we re-write: 

5B _ I 2B 

m" iF If- 

We can face this expression as a first order differential equation whose most general solution is given 

A = -eS^. (10) 

where C is a constant. From this expression we find the following algebraic equation: 

n2- + n-C7 = o, (11) 



where C is a constant. The solutions are: 



^ 2A ' 

The constant C is chosen in such a way that, in the hmit of symmetric nuclear matter, the polarizability 
n is the one calculated in expression (|5|). This leads to C = —p/{4:Asym)- 

A safer and sounder way of obtaining solutions for the polarizability is to find the value of [3 at 
which expression (P) is satisfied. We obtain a quadratic equation for (3 whose solutions are given by: 



' 5(3 V ^(0,1)2 5[3 

There are two solutions for this expression and we will be concerned only with one of them (which 
produces sounder results) in the limit of very small amplitude e. In this limit for the positive sign of 
expression (|l^) the solution is trivially found to be: 

n.=o,t=i = ^ ^ — p (14) 

This expression has exactly the same form of symmetric nuclear matter. The difference is that, here, 
A does depend on b. 

It is also possible to extract an explicit isospin dependence of the generalized symmetry energy 



coefficient. From the solution of the polarizability (14) we consider the first derivative with relation 
to 6: 

d/3 ep dA f3 dA ,^ ^.x 

lib ^ Z^lh ^ ~Alh' ^ ^ 

A complementary expression can be obtained from the relation between b and /3 of (^. It yields: 

dp p 



db 2 + b' 
Equating these two last equations we obtain: 



(16) 



From which it is possible to derive the following relation between the isospin s.e.c. and the n-p 
asymmetry: 

A = Ay^(^^±^]\ (18) 
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where Z is a constant which depends on the channel (isovector in the present case) and eventually on 
other properties of the system and Asym is the s.e.c. of symmetric nuclear matter. In the symmetric 



case 6 = and hence A = Agym- For 6 = 2 (neutron density three times larger than the proton 
density) we obtain A = 2^ Asym- If we consider Z = 1 a considerably high value is obtained. It is 
therefore missing some dynamical information to fix the constant Z for each channel. In the next 
section we find that different Skyrme forces will result different values of Z. We want to emphasize 
that the prescription (Q) was the relevant information for this calculation. Any other relation between 
b and /? will induce different asymmetry dependence of the symmetry energy coefficients. 

Besides that expression ( [l^ corresponds to the static limit of the dynamical polarizability [17|. It 
defines the static screening functions A^'^ in the usual way: 

e = C'A^'^P, 

where C = 1/po- From this expression we see that the external field e induces density fluctuations 
P according to the intrinsic properties of the medium, given by the screening function C'A^'^. These 
arguments are valid for any perturbation of the other channels for asymmetric nuclear matter in isospin 
as well as in scalar channels yielding the functions A^'*^. 

We have shown that the polarizabilities of asymmetric nuclear matter yield the generalized symme- 
try energy coefficients. The dependence of the s.e.c. on the n-p asymmetry can be deduced unless for 
a constant exponent Z. It is however necessary to provide a prescription for the induced fiuctuation 
p. A nearly exact expression for the dynamical polarizability of a non relativistic hot asymmetric 
nuclear matter for variable densities was derived with Skyrme interactions in [|2^. Its static limit is 
exploited in the following sections for more general Skyrme forces. 

2.2 Screening functions with Skyrme forces 

In 17| the dynamical polarizabilities of asymmetric nuclear matter were calculated using Skyrme 



forces. An alternative and independent calculation for the dynamical polarizability was done in [l24| . 
Although the static limit has been found yielding the screening functions it has not been exploited 
sufficiently well. It can also be derived directly from the static Hartree-Fock equation with the same 
external perturbation. The static generalized symmetry energy coefficients of hot asymmetric nuclear 
matter, As,t, can be written as: 

As,t = -7^4 1 1 + '^WNc + 6F/'*m;(/5, + pa) + 12M;v,'^'vf {N,pa - PcN^) + 
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This expression was derived in ||17|] and is one of the main concerns of the present article The 
densities p^, and are given by: 

Pv = vpn + (1 - v)pp, 

= vMn + (1 - v)Mp, 
Ny = vNn + (1 - v)Np, 

where v stands for n-p asymmetry coefficients (c, d) defined below (a measure of the fraction of neutron 
density). The above densities are defined by: 

2M* r 

iNg,Pg,M,) = / dfgik)ik,k^k^). 



In these expressions dfq{k) are the differential fermion occupation numbers for neutrons {q = n) and 
protons {q = p)- In the zero temperature limit they reduce to delta functions at the Fermi surface 
momentum dfq — > —6{k — kp). At non zero temperature one has to fix the chemical potential to 
calculate the Fermi occupation number and the densities. But this study of the s.e.c. dependence on 
the temperature would be more appropriately done with an equation of state and it is beyond the 
scope of this article. Vq and Vi are functions of the Skyrme forces parameters (calculated for Skyrme 
forces with more general density dependence and shown in the appendix) and M* = m*/(l + a/2) is 
an effective mass for the proton. Besides that, the four asymmetry coefficients are: 

ml pop 2 + 6 1 + (1 + 5)3 

(The coefficient b is related to a frequently used asymmetry coefficient: a = (2pon — Po)/ Po-, by the 
expression: 6 = 2a/(l — q)). 

In the scalar channel the screening function will be referred to as the dipolar incompressibility, 
K£), and is related to the usual nuclear matter incompressibility, K^a by the following expression in 
terms of the SLyb Skyrme force parameters 



KD = K^ + ^TF-2ViklpQ + ^Up'^+^. (21) 



This relationship is slightly different for the SkSC interactions In this reference the authors fitted 
neutron star properties with parameterizations that yielded different values for the n-p symmetry 
energy at saturation density: SkSC4 {27MeV), SkSC6 (30MeF) and SkSClO (32Mey). At different 
densities and n-p asymmetries (and also temperatures) the static screening functions have different 
values. They are studied in the next section. 
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3 Results and discussion 



In this section we show some figures which exhibit the isospin and density dependences of the nuclear 
matter screening functions of expression (|^) for the isospin and scalar channels. For this we use the 
functions Vq and Vi (written in the appendix) calculated for Skyrme functions with different density 
dependences. The parameters of one of them (SLyb) were fitted from results of neutron matter 
properties obtained from microscopic calculations in [^]. Other forces (SkSC4, SkSC6 and SkSClO), 
which have slightly different density dependences, had their parameters fixed by adjusting a large 
amount of nuclear masses yielding the same results of the Extended Thomas Fermi method and shell 
corrections calculated by the Strutinsky-integral method ||l^, . One of the strong characteristic of 
the forces SkSC is the value of the effective mass which is kept to be equal to the free nucleon mass. 
This may be interpreted as a result of the coupling between particle modes and surface vibration 
modes. 

In figure 1 the generalized isospin symmetry energy coefficient ^0,1 is shown as a function of the 
ratio of the density to the saturation density for Skyrme interactions SLyb, SkSC4, SkSC6 and SkSClO. 
Since Skyrme forces are not necessarily expected to describe physics at high densities we decided to 
investigate the dependence of the screening functions up to 1.5 po- Most of the isovector screening 
functions of figure 1 does not vary much. Most of them have the tendency to keep nearly constant in a 
value close to that at the saturation density and decrease at higher densities. For the SLyb force two 
curves are shown. ^0,1 corresponds to the complete s.e.c. expression (19) whereas Ur neglects the Vi 



term present in expression ( A.l ). They exhibit nearly the same behaviour. The SkSClO force, however, 
yields a different behaviour which is more compatible with that of the parameterizations of expression 
(H), with F{u) = u from (plotted with circles). There are no means (not yet) of choosing the most 
realistic (and reliable) of the forces in this channel for varying densities. The force SkSClO reproduces 
rather the behaviour of relativistic parameterizations. In figure 2 the generalized isospin s.e.c. ^0,1 
is shown as a function of the n-p asymmetry b for the different forces at the saturation density po- 
There is a trend which is absolutely dominant for all the interactions. It is the increase of .Ao,i with 
b. One big difference between these curves for each of the forces is given by the slope which can be 
larger (eg. for SkSC6) or smaller (for SLyb). It seems to us that if the slope were so large as that 
obtained from SkSC6 it would have already manifested in studies of neutron matter as well as in heavy 
nuclei observables (for the lead nucleus b = 0.54). Further studies of neutron stars, supernovae and 
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heavy ions collisions may lead to a more precise determination of this slope. In section 2.1 a general 



expression (18) for the isospin dependence of the isospin s.e.c. has been derived making use of relation 
d^) . However it was not possible to fit any of the results obtained using the above Skyrme forces with 
expression (^), for any Z. The form of the curves (slopes) obtained with the used Skyrme forces, in 
figure 1, indicates that Z would be smaller than one whereas the numerical values of ^o,i (with the 
same Skyrme interactions) suggest that Z > 1. This indicates that the calculation done in section 2 
is not completely compatible with the Skyrme forces which have been used in this paper. This could 
eventually be solved by means of microscopical calculations (e.g. using Monte Carlo techniques as 
done in [^]) and could, hopefully, also be checked by new parameterizations of mass formulae. 

Another issue which can be taken from figures 1 and 2 is the isospin dependence of the s.e.c. at 
diverse densities. The effects of changing the asymmetry and the density sum up (with relation to 
the symmetric nuclear matter at saturation density) reproducing behaviours already present in figures 
1 and 2 for given p and b. This means that the same results of figure 1 (density dependence) will 
be present for nuclear matter with an asymmetry b (shown in 2). The difference is that: the more 
asymmetric is the nuclear matter more repulsive is the isovector interaction (and hence higher the 

A,i)- 

In ||l9| the EOS of a collapsing star is studied and the results support forces for which is 
considerably greater than 27MeV to keep neutron matter stable. However these higher values spoil 
the fit of nuclear masses within the frame of the extended Thomas- Fermi plus Strutinsky integral mass 
formula Q. Our results make these discrepancies be compatible in a certain sense. We have found 
that for larger values of the n-p asymmetry b the isospin s.e.c. increases considerably corresponding 
to the case of neutron stars medium. Therefore for the same Skyrme force several values of Ur, or 
rather Aq^i = Ao^i{b,p), can be considered. For the description of (not very asymmetric) nuclei the 
lower values can be expected to yield better results. 

For the generalized scalar dipolar s.e.c, figures 3 and 4 show that, roughly speaking, all the 
interactions are attractive for lower densities (figure 3) and higher n-p asymmetries (figure 4). This 



corresponds to instabilities which are related to those found, for example, in |26]. For neutron stars 
it may seem that all these forces (except SLyb) lead to instabilities according to figure 4 (for fixed 
density po). However higher densities may stabilize the system (while lower densities lead to still more 
unstable systems) . As discussed for the isospin case the effects of varying the nuclear density and the 
asymmetry can be considered simultaneously, i.e., to obtain the n-p asymmetry dependence of the 
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scalar s.e.c. at a density p Pq the behaviours of figures 3 and 4 can be analysed together. Therefore, 
for higher densities (making the ^0,0 to increase) the system may experience higher asymmetries 
(6) without the instabilities present for forces SkSC4 and SkSC6. We do not exhibit, however, a 
quantitative estimate of the instabilities points in the variables h and p due to the diverse behaviour of 
Skyrme forces. Moreover these points may be out of the validity of the application of these interactions. 
For the moment there is no way of verifying this. The interaction in this channel becomes more and 
more repulsive for higher densities. The increase of the n-p asymmetry on the other hand, make the 
interaction continuously more attractive for most of the forces as shown in figure 4. However, the 
slopes for interactions SkSC seem to be too negative. In figure 4 it is also shown the parameterization 



for n-p asymmetry dependence of the incompressibility of nuclear matter {Kao{h)) of reference [27 
-from the equation of state- (not the incompressibility K^o itself, but its asymmetry dependence- for 
the case of SLyb). This parameterization (using the n-p asymmetry parameter of our work) is given 
by: 

K{h) = K{b = Q){l-a'{b/{h + 2)f) (22) 

where a' (in ^^) is of the order of 1.28 up to 1.99 for several Skyrme interactions (zero temperature) 
and K{b) stands for Kao or ^0,0 = K^. It is seen in figure 4 that this case (from [^) has nearly 
the same behaviour (dependence on the n-p asymmetry b with p = po) as our calculation for force 



SLyb, a remarkable feature. In another work |19| it was found the same kind of dependence on the 



n-p asymmetry for the incompressibility of nuclear matter but with other numerical coefficients: 

K{h) = K, + Ks{h/{h + 2)f (23) 

These parameters assume the following values for the forces used in that article: 

SkSCA{K^ = 234.7MeV, = -334.9MeF); SkSCl^{K^ = 235AMeV, = -203.5Mey); 
SkSCQ{K^ = 235.8MeV, Ks = -136.6MeF). 

(24) 

These values are (a little smaller) of the same order of those obtained in the figure 4 and expression 



(22). It can expected that this dipolar incompressibility is also of relevance for the study of the 
Isoscalar Dipole Giant Resonance ]2^ ]. 

Consequences concerning the supernovae mechanism can be discussed now. The higher the n-p 
symmetry energy coefficient the smaller is the deleptonization (electron capture) in the quasi-static 
phase of the supernovae mechanism yielding a larger final proton fraction and faster cooling (via 
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neutrino emission). As we have found that ^0,1 increases with the n-p asymmetry it can be expected 
that as deleptonization proceeds it becomes more and more difficult to convert protons into neutrons. 
This also means that the final (neutronized) star would be less asymmetric than expected with smaller 
(fixed) values for Ur- This picture yields a stronger shock wave since the energy loss inside the star 
due to deleptonization is smaller. Therefore the increase of isospin symmetry energy coefficient (which 
occurs for higher n-p asymmetries) helps a successful explosion of the (contracting) star (keeping the 
density fixed). This conclusion has the essential feature compatible with the analysis done for the SN 
1987a event g, |9|,|3|. 



4 Conclusions 

Summarizing, generalized symmetry energy coefficients of nuclear matter were investigated. Firstly 
we have shown that the polarizabilities of asymmetric nuclear matter yield the generalized symmetry 
energy coefficients. The dependence of the isospin s.e.c. on the n-p asymmetry, by means of the 
relation (|^) was found to be: 

"'^. c 



^sym 



[2 + by ^ 
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where Z is a constant. To obtain such expression it was necessary to provide a prescription for the 
dependence of the induced fluctuation f3 on the n-p asymmetry b (expression (7)) p3[| which could be 
chosen differently. The form of the curves obtained from the used Skyrme forces in figure 1 indicate 
that Z would be smaller than one whereas the numerical values of ^0,1 suggest that Z > 1. This 
constant seems to contain relevant dynamical information and it could be determined by the study of 
the symmetry energy coefficient microscopically. 

Secondly a study of the generalized symmetry energy coefficients (isovector and scalar) was done 
with Skyrme forces. However these results are not adjusted by the function obtained (the expression 



(18)) The use of Skyrme-type interactions allowed to obtain analytical expressions for the s.e.c. Their 
density and n-p asymmetry dependences were analyzed for different Skyrme forces which may yield 
very different behaviours. In some cases it is possible to discard unphysical results but there are 
results which cannot be chosen due to non existing experimental knowledge. Nevertheless neutron 
stars and eventually very asymmetric nuclei properties as well as heavy ion collisions can provide 
valuable information. Although one should not believe that only one Skyrme force parameterization 
could account for the description of all nuclear observables at different n-p asymmetries (as well as 
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densities and temperatures) it is acceptable the idea that several parameterizations could hopefully 
describe different ranges of the dependence of nuclear observables (as the s.e.c.) with these three 
variables. Nevertheless it is tempting to conclude that the force SLyb, which was adjusted to reproduce 
general properties of asymmetric nuclear matter, yielded, in general, more reasonable behaviour for 
the symmetry energy coefficients for varying n-p asymmetries. The dependence of the isospin s.e.c. 
on the n-p asymmetry indicates that the deleptonization process in supernovae is suppressed as the 
neutronization occurs. This would make final proton fraction bigger with less energy loss helping the 
SN explosion (and also less asymmetric neutron stars matter) in agreement with expectations. A 
quantitative estimate of this effect however will only be possible after a more precise determination of 
the dependence of the s.e.c. on b. 
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Appendix: Functions F/'^ 

We write the parameters Vq and Vi from the linearization of the Hartree Fock equation (for nuclear 
matter) for the Skyrme forces SkSC used by [19|. They contain more involved density dependence: 

1 - ' m*po^p 

where only the Vq and V2 were modified. 

In the SkSC forces the density dependence implies that this part of the interaction between two 
protons depends only on the proton density being more reasonable. It can be written in the following 
form: 

|(1 + xM {ai{p,, + p,jrS{r,,j) + asp") , (A.2) 

where is the two-body spin exchange operator, q indicates proton or neutron and oi, 02 are taken 
to be 0, 1 or 1,0 for Skyrme interactions SkSC and SLyb respectively. 
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In the isovcctor channel one obtains: 



0,1 



to 
2 



^0 + ^1- 



12 



a2 (x3 + ^) + «i (l + y) - - ^3)(a + 2)(a + 1) 



^1°'' = i7:{t2{l + 1x2) - ti(l + 2xi)), 



0,1 



16 

h 



a2(2+2^3)aPo ^(c/9„ + (c - l)pp) + 

X3, 



+ai (1 + (cpn + Pp{c - 1)) + 2(1 - x^){a + 2)(a + l)(cp^ + p^(c - 1))/16 



/12, 



(A.3) 



where poru Pop and are the proton, neutron and total saturation densities of asymmetric nuclear 
matter. For the scalar channel: 



^0,0 



0,0 



(^3^ + (« + l)(a + 2)t3P^ 
3^ + (5 + 4..4 



ai(l + f)( 



2T6j 16+"^(' + T)i2 



(l + 6c), 



0,0 



^3 

12 



{xz + .^){cpn + {c-1)PpPq ^'^) + 



+ai(a + l)(a + 2)(l -X3 



(2 + 6)"+2 ^ ^ (2 + 6)°+2 



/2-a2(l + (l + 6)V)/(2 + 6) 



(A.4) 
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Figure caption 



Figure 1 Isovector screening function Aq^i = p/{2Ilj^ ) of symmetric nuclear matter as a function 
of the ratio of density to density at saturation (n = p/ po) for interactions SLyb (solid), SkSC4 (dotted), 
SkSC6 (dashed), SkSClO (long dashed). The simplified expression (a,-), i.e. without terms of order of 
Vi, for the force SLyb (thin dotted-dashed). Circles (P.A.L. 1) for the expression of ar as a function 
of u with three different functions F{u) from reference p^ . 

Figure 2 Isovector screening function Aq^i = p/(2n^^) as a function of the asymmetry coefficient 
b (at pq): dotted line for SkSC4, dashed lines for SkSC6, solid line for SLyb and long-dashed lines for 
SkSClO. 

Figure 3 Scalar (dipole) screening function Aq q = p/{2Il^jf) as a function of the ratio of density 
to density at saturation (u = p/po) for the different interactions with the conventions of figure 1. 

Figure 4 Scalar (dipole) screening function Aq^ = p/{2Il^^) as a function of the asymmetry 
coefficient b, with pQ, with the different forces (the conventions of figure 2) and also the n-p asymmetry 
dependence of the compressibility of nuclear matter of reference using K£){b = 0) from the 
SLy force (star). 
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